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ABSTRACT 

NAVSYS  has  developed  a  GPS  receiver,  the  HAGR  (High- gain  Advanced  GPS  Receiver), 
which  uses  digital  beam  steering  to  increase  the  precision  of  the  GPS  observations.  In  this  paper, 
the  HAGR  is  described  and  test  data  is  included  showing  its  enhanced  performance  for  high 
accuracy  scientific  applications. 

When  operating  in  a  differential  or  kinematic  positioning  mode,  HAGR  provides  low  noise 
observations  of  the  pseudo- range  and  carrier- phase  data.  By  coherently  combining  signals  from 
as  many  as  16  antenna  elements,  the  gain  to  each  satellite  tracked  can  be  increased  by  over  10 
dB.  Moreover,  the  digital  beam  steering  approach  used  also  acts  to  reduce  the  errors  contributed 
by  multipath  effects.  Test  data  is  presented  on  the  HAGR  measurement  accuracy  for  high 
accuracy  scientific  applications. 

In  this  paper,  the  GPS  carrier- phase  time  transfer  technique  is  described  and  a  discussion  is 
included  on  the  error  components  that  currently  limit  the  time-transfer  accuracy  using  this 
method.  Sub -nanosecond  time  transfer  capabilities  are  expected  using  GPS  carrier-phase 
observations.  Testing  to  date  has  shown  that  conventional  GPS  receivers  introduce  significant 
time  offsets  in  the  carrier  phase.  These  offsets  are  on  the  order  of  1-2  nanoseconds  and  currently 
dominate  the  error  budget  when  performing  carrier-phase  time  transfer.  Test  data  is  also 
presented  from  the  HAGR  showing  the  accuracy  of  the  code  and  carrier  phase  observations  for 
time  transfer  applications. 
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INTRODUCTION 


The  HAGR  design  is  based  on  NAVSYS’  Advanced  GPS  Receiver  (AGR)  PC- based  digital 
receiver  architecture  integrated  with  a  digital  beam  steering  array,  see  Figure  1[1] .  Using  a 
proprietary  digital  signal  processing  algorithm,  the  HAGR  is  able  to  combine  the  GPS  signals 
from  as  many  as  16  antennas  and  create  a  multi- beam  antenna  pattern  to  apply  gain  to  up  to  eight 
GPS  satellites  simultaneously.  The  HAGR  16-element  antenna  array  is  shown  in  Figure  2. 

Testing  has  also  been  conducted  with  our  miniature  antenna  array,  shown  in  Figure  3  [2,3]. 


Figure  2  HAGR  16-element  antenna 


Figure  1  High-gain  Advanced  GPS 
Receiver  (HAGR) 


array 

The  performance  specifications  for  the  HAGR  for  a  16- element,  LI  C/A  code  version  of  this 
product  are  included  in  reference^.  Currently  an  L1/L2  Precise  Position  System  (PPS)  version 
of  the  HAGR  (the  HAGR- 200)  is  also  in  development. 


A 

Figure  3  6”  diameter  4-element  Mini-Array 


HAGR  SYSTEM  ARCHITECTURE 


16  Antenna  Elements  10  Processin8  Channels 


Figure  4  HAGR  System  Block  Diagram 


The  HAGR  system  architecture  is  shown  in  Figure  4.  The  signal  from  each  antenna  element  is 
digitized  using  a  Digital  Front-End  (DFE).  The  bank  of  digital  signals  is  then  processed  by  the 
HAGR  digital- beam- steering  card  to  create  a  composite  digital  beam- steered  signal  input  for 
each  of  the  receiver  channels.  Figure  5  illustrates  the  beam- forming  satellite  geometry  and 
mathematical  expression  for  line  of  sight  (los). 
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Figure  6  16-element  array  composite 
beam  pattern 


Figure  5  Beam  forming  satellite  geometry 


If  attitude  data  (pitch,  roll,  yaw)  is  provided  from  an  inertial  navigation  system  or  attitude  sensor, 
the  HAGR  will  operate  while  the  antenna  is  in  motion[5].  The  default  mode,  for  static  operation, 
is  to  align  the  array  pointing  north. 

The  digital  beam  forming  provides  significant  benefits  in  improving  the  measurement  accuracy 
due  to  the  narrow  beam  antenna  pattern  directed  at  each  satellite  tracked.  As  shown  in  Figure  6, 
a  16- element  array  will  provide  up  to  12  dB  of  additional  gain  on  each  satellite  tracked. 


DGPS  AND  KGPS  NAVIGATION  ACCURACY 


The  accuracy  of  a  differential  GPS  (DGPS)  solution  is  a  function  of  the  solution  geometry  and 
the  accuracy  of  the  raw  pseudo- range  measurements  as  shown  in  Equation  1  and  Equation  2.  A 
similar  equation  also  exists  for  a  Kinematic  GPS  (KGPS)  solution  which  is  derived  from  carrier 
phase  observations  and  a  solution  of  the  integer  carrier  cycle  ambiguities  (see  Equation  3). 
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Equation  3 


E[xTlGPS xtGPS  ]  =  PDOP 2a  2gps  *  PDOP 2 (a  2cph  +  a 2Mcph)l2 

The  major  benefits  of  the  HAGR  digital  beam  forming  for  DGPS  and  KGPS  navigation 
applications  are  the  increase  in  accuracy  of  the  pseudo- range  and  carrier-phase  observations  and 
the  reduction  in  multipath  errors  due  to  the  antenna  directivity.  These  benefits  are  discussed  in 
this  section. 

The  differentially  corrected  pseudo- range  accuracy  is  dominated  by  two  error  sources  as 
illustrated  in  Equation  4.  The  first  is  receiver  noise  and  the  second  is  the  multipath  error.  The 
receiver  noise  is  a  function  of  the  effective  delay- lock- loop  bandwidth  after  carrier  smoothing  is 
applied.  This  can  be  computed  from  the  following  equation  where  Tc  is  the  C/A  code  chip 
length  (293  meters),  d  is  the  correlator  chip  spacing  and  C/No  is  the  received  signal- to- noise  ratio 
in  dB-Hz. 

Equation  4 
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This  is  plotted  in  Figure  7  against  C/No  assuming  a  1  Hz  bandwidth  delayed  locked  loop  (no 
carrier  smoothing),  a  code  chip  separation  of  1  and  an  accumulation  frequency  of  1-kHz. 


Figure  7  Sigma  PR  Noise  vs  C/NO  (Bdll=l,  d=l) 


From  Figure  7,  a  10  dB  increase  in  signal  power  will  have  the  affect  of  reducing  the  pseudo¬ 
range  noise  by  a  factor  of  0.3.  The  HAGR  routinely  receives  satellite  signals  with  C/No  of  54  to 
58  dB-Hz  which  reduces  the  pseudo- range  noise  level  to  below  10  cm  with  carrier  smoothing. 

HAGR  MEASUREMENT  ACCURACY 

The  increase  of  about  10  dB  in  signal- to- noise  ratio  from  the  HAGR  digital  beam  forming  is 
shown  in  Figure  8  to  Figure  10,  compared  against  two  conventional  GPS  reference  receivers'1  ]. 
From  these  plots,  it  can  be  seen  that  the  HAGR  C/No  is  significantly  higher  than  the  reference 
receiver,  demonstrating  the  effect  of  the  gain  from  the  digital  beam  forming. 
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Figure  8  SNR  Comparison  Between  16- 
Antenna  HAGR  and  Novatel’s  for  PRN  2 
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Figure  9  SNR  Comparison  Between  16- 
Antenna  HAGR  and  Novatel’s  for  PRN  3 
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Figure  10  SNR  Comparison  Between  16-Antenna  HAGR  and  Novatel’s  for  PRN  13 

The  increased  gain  also  results  in  improved  pseudo- range  and  carrier- phase  tracking  performance, 
and  the  directionality  of  the  beam- steering  antenna  array  reduces  the  effect  of  multipath  on  the 
solution.  In  Table  1,  the  short  term  noise  is  listed  for  each  of  the  two  HAGR  units  tested.  The 
gain  provided  by  the  beam  steering  has  maintained  the  signal- to- noise  generally  above  50  dB-Hz, 
providing  sub-meter  level  short  term  noise  on  the  pseudo- range  performance.  This  increased 
accuracy  reduces  the  time  needed  to  resolve  the  carrier- cycle  ambiguities  needed  for  computing 
kinematic  position  solutions  or  a  carrier-phase  time  transfer  solution. 

Table  1  HAGR  PR  Noise  Performance  Data 
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Multipath  errors  are  caused  by  the  receiver  tracking  a  composite  of  the  direct  GPS  signals  and 
reflected  GPS  signals  from  nearby  objects.  The  resulting  pseudo- range  error  and  carrier-phase 
error  is  a  function  of  the  phase  offset  between  the  direct  and  multipath  signals  and  the  relative 
signal  strength.  For  a  fixed  installation,  these  errors  appear  as  biases,  changing  only  as  the  line- 
of- sight  to  the  satellite  changes  due  to  the  satellite  motion.  Multipath  mitigation  techniques  have 
been  developed  using  multi- correlator  techniques  that  improve  the  performance  of  the  code 
tracking  loops  in  the  presence  of  multipath.  However,  these  have  little  effect  against  carrier- 
phase  errors  introduced  by  the  multipath  signals.  The  HAGR  digital  beam- former  has  the 
advantage  that  the  multipath  is  reduced  on  both  the  pseudo- range  and  the  carrier-phase  errors 
through  the  directivity  of  the  antenna  pattern  towards  the  satellite  which  reduces  the  effect  of  the 
multipath  signal. 


The  effect  of  multipath  on  the  GPS  signals  can  be  modeled  through  the  following  equations. 
Equation  5 

s(?)  =  AC(t  +  x)  sin(cor  +  0 )  +  A^C(t  +  x  ^ )  sin((D?  +  0^ ) 
s(t)s(t)  =  AR(t  -f)sin(cor  +0)+  AmR(tm  -f)sin(cof +  0m) 

=  (y4  +  A)sin(co?+0  +0) 

The  above  equation  can  be  solved  for  the  pseudo- range  error  that  will  be  observed  with  the  DLL 
tracking  loops  and  the  phase  error  that  will  be  observed  by  the  PLL  tracking  loops.  This 
simplifies  to  the  following  expression,  if  it  is  assumed  that  the  multipath  reflections  are  relatively 
close  to  the  antenna  (compared  with  the  C/A  code  chip  length  of  293  m)  and  that  R(t)  from  the 
DLL  is  approximately  equal  to  one  for  both  the  direct  and  multipath  signals. 

Equation  6 


This  equation  represents  the  approximate  maximum  error  for  a  1-chip  Early/Late  correlator 
assuming  that  the  multipath  delays  are  less  than  Vi  chip.  If  the  multipath  errors  are  10  dB  down 
from  the  satellite  signal  power,  then  the  multipath  error  will  be  approximately  10%  of  the 
multipath  delay.  For  example,  if  signals  are  received  with  a  delay  of  20  meters,  then  the 
multipath  error  will  be  less  than  2  meters  on  the  pseudo- range  observation.  Since  the  beam¬ 
forming  antenna  provides  10  dB  gain  in  the  direction  of  the  satellite  signal,  the  multipath  signals 
are  received  with  at  least  10  dB  lower  power  than  the  direct  satellite  signals. 

Equation  7 

AIA=\A  +  AMe^\lA 

0~  =Z(A  +  AMe* ) 

A0  =0-0M 

In  Figure  1 1  the  effect  of  the  multipath  signal  on  the  C/N0  envelope  is  shown  as  a  function  of  the 
multipath  signal  power.  In  Figure  12  the  multipath  phase  angle  errors  and  signal  amplitude  are 
shown  versus  the  multipath  phase  angle  offset  (A0)  for  the  cases  when  Am=A  (0  dB),  Am=A/V2 
(-3dB  down) ,  Am=A/2  (-6dB  down),  and  Am=A/V  10  (- 10  dB  down).  This  figure  illustrates  the 
benefit  of  the  beam  former  in  also  limiting  multipath  errors  on  the  phase  observations  -  a  key 
area  of  concern  for  kinematic  GPS  applications. 


Figure  11  Multipath  Amplitude  Effect  Figure  12  Multipath  Phase  Angle  Error 


The  multipath  effect  on  the  GPS  measurements  is  difficult  to  observe  directly.  The  best 
indication  of  the  magnitude  of  multipath  errors  is  from  the  variation  in  C/NO  and  the  short-term 
difference  between  PR+CPH  which  are  both  dominated  by  multipath  errors.  Test  data  taken 
with  two  HAGRs  indicated  that,  with  beam- forming,  the  majority  of  the  satellite  signals  have  an 
C/NO  delta  of  less  than  2  dB  [4],  From  Figure  1 1  and  Figure  12,  the  carrier  phase  variation 
would  be  expected  to  be  less  than  +/-  5  mm  when  the  multipath  signal  is  below  this  level. 

It  should  be  noted  from  Figure  12  that  a  strong  multipath  signal  can  cause  phase  variations  on  the 
order  of  +/-  2  cm  which  would  severely  affect  the  kinematic  GPS  performance  for  a  single¬ 
element  GPS  receiver.  To  compare  the  performance  of  the  antenna  anay  with  and  without  the 
beam  steering,  the  residual  phase  enor  on  each  of  the  individual  elements  of  the  antenna  array  is 
plotted  in  Figure  13.  This  shows  that  the  carrier  phase  error  can  easily  be  offset  by  +/-  0.1  cycles 
(2  cm)  between  the  individual  elements.  For  high  accuracy  applications,  the  multipath  rejection 
from  spatial  processing  the  HAGR  array  signals  can  significantly  reduce  this  error  source. 
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Figure  13  Carrier-phase  multipath  without  beam  steering 
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Figure  14  PR+CPH  -SV  21 


Figure  15  PR+CPH  -SV  23 


KINEMATIC  GPS  TEST  RESULTS 

The  kinematic  performance  of  the  HAGR  antennas  was  tested  by  setting  each  of  the  antennas  on 
two  survey  marks.  The  remote  antenna  array  was  located  0.793  meters  East  of  the  reference 
antenna  array.  The  NAVSYS’s  kinematic  GPS  software  was  used  to  process  the  data.  Figure  16 
to  Figure  1 8  show  the  processing  results.  During  the  test,  6  valid  satellites  were  available.  These 
test  results  show  that  the  kinematic  GPS  positioning  error  achieved  a  standard  deviation  of  4  mm 
(1-sigma)  in  the  north  and  east  directions  and  6  mm  (1- sigma)  vertically.  This  is  consistent  with 
a  carrier  phase  measurement  accuracy  of  3  mm  (1-sigma).  These  test  results  show  that  the 
multipath  errors  on  the  carrier  phase  are  maintained  on  the  order  of  a  few  millimeters  by  the 
HAGR  beam  forming.  The  mean  offset  between  the  two  antennas  observed  by  the  kinematic 
solution  was  0.7933  meters  which  agreed  to  the  surveyed  location  to  within  the  precision  of  our 
measurement  capability. 


Figure  16  HAGR  DGPS  Position  Solution 


Figure  17  HAGR  KGPS  Position  Solution 


Figure  18  HAGR  Kinematic  Position  Solution  (NED) 

GPS  CARRIER  PHASE  TIME  TRANSFER 

In  this  section  of  the  paper,  the  performance  of  the  HAGR  in  performing  GPS  carrier-phase  time 
transfer  technique  is  described.  Previous  testing  with  conventional  GPS  receivers  has  shown  that 
carrier  phase  instabilities  can  cause  offsets  on  the  order  of  1-2  nanoseconds,  which  currently 
dominates  the  error  budget  when  performing  carrier- phase  time  transfer.  The  ability  of 
NAVSYS’  HAGR  to  make  phase  coherent  measurements  from  multiple  antenna  elements,  also 
means  that  it  provides  a  highly  phase  stable  observation  of  the  GPS  carrier,  relative  to  a  local 
reference  oscillator.  This  allows  the  HAGR  to  provide  highly  precise  time  observations  for 
carrier- phase  time  transfer. 

GPS  carrier- phase  measurements  provide  the  potential  for  much  improved  precision  in  time  and 
frequency  transfer  [6,7,8].  Time- Transfer  errors  approaching  100  picosecond  (ps)  are  expected 
using  this  approach.  Typical  carrier  phase  measurement  noise  can  be  on  the  order  of  ten  ps 
whereas  the  code  measurement  noise  can  be  as  high  as  ten  nanoseconds  (ns).  Multipath  errors 
are  also  much  smaller  on  the  carrier-phase  observations  than  on  the  code-based  pseudo-range 
measurements.  Accordingly,  GPS  carrier-phase  measurement  accuracy  is  100  to  1000  times 
better  than  the  code  based  pseudo-range  measurements. 

The  GPS  system  errors  that  affect  the  accuracy  of  the  carrier-phase  time  transfer  performance  are 
listed  in  Table  2 


Table  2  GPS  System  Errors 

1.  Dual  Frequency  Ionosphere  errors  (calibration  bias,  increase  noise) 

2.  Troposphere  errors  (Weather  data,  Models) 

3.  Receiver  Measurement  Noise 

4.  Multi- path  Noise 

5.  Satellite  Position  (Orbits) 

6.  Station  Position  (Location) 


Test  results  have  shown  that  the  dominant  errors  currentiy  affecting  the  accuracy  of  carrier-phase 
time  transfer,  are  not  the  GPS  system  errors  shown  in  Table  2,  but  are  due  to  environmental 
effects  within  the  GPS  receiver. 


The  time  delay  of  the  GPS  signal  as  it  propagates  through  a  complete  GPS  receiving  system 
consists  of  the  delay  through  the  GPS  receiver,  GPS  antenna  cables  and  the  GPS  antenna  with  its 
associated  antenna  electronics.  All  of  these  GPS  receiving  system  sub-components  are  affected 
by  environmental  influences.  Studies  of  the  temperature  sensitivities  of  several  of  these  GPS 

receiving  systems  have  shown  delay  variations  of  as  great  as  several  nanoseconds  per  degree  C.[9 
10  11] 


Table  3  GPS  Receiver  Temperature  Sensitivity 


Temperature  effect 

Receiver  Code 

Measurements 

(150-  1500)  ps 
per  °C 

Receiver  Carrier 
Measurements 

(10  -  200)  ps  per 
°C 

Antenna  cable 

0.5  ps  per  °C  per 
Meter 

Antenna  electronics 

(5  -  50)  ps  per  °C 

Since  all  of  these  temperature  effects  are  common  to  all  receiver  channels,  these  errors  are 
mapped  into  the  users  local  clock  error.  This  does  not  affect  the  use  of  this  data  for  typical  geo¬ 
location  application,  but  for  time  transfer  applications  these  temperature  effects  must  be 
minimized.  Specially  constructed  phase- stabilized  antenna  cables  can  be  used  that  will  reduce 
the  delay  fluctuations  through  the  antenna  cable  by  a  factor  of  20  or  more.  However,  the  GPS 
receiver  front-end  itself  must  also  be  designed  to  provide  a  highly  stable  carrier-phase  reference 
over  temperature  variations.  This  has  been  accomplished  in  the  HAGR  design. 

TIME  TRANSFER  LAB  TEST  RESULTS 

To  test  the  time  transfer  performance  of  the  HAGR  receiver,  two  receivers  were  set  up  to  operate 
using  a  common  10  MHz  time  reference  and  also  a  common  antenna.  This  test  will  cancel  the 
GPS  system  errors  shown  in  Table  2,  leaving  the  effect  of  the  carrier  phase  observation  and 
uncalibrated  receiver  errors  on  the  solution. 

The  raw  carrier  phase  difference  was  computed  between  the  two  receivers  for  each  satellite 
tracked.  This  was  corrected  for  the  integer  ambiguity  offset  only.  The  residual  error  between 
two  data  sets  for  each  satellite  is  plotted  in  Figure  19  and  Figure  20.  The  HAGR  was  power- 
cycled  between  these  two  data  sets.  As  can  be  seen,  both  data  sets  observed  a  common  bias 
between  the  units  of  around  0.02  cycles  and  has  a  standard  deviation  of  the  carrier- phase 
difference  residual  of  16  psecs.  Each  satellite  observes  a  common  offset  between  the  units  of  14 
psecs  +/-  3  psecs,  indicating  that  the  HAGR  units  should  be  able  to  be  calibrated  to  this  level  by 
averaging  the  satellite  observations. 


Table  4  Carrier-phase  time  difference  accuracy 


SVID 

1 

14 

16 

18 

22 

25 

Mean  offset 
(cycles) 

0.022 

0.020 

0.022 

0.026 

0.022 

0.020 

Mean  offset 
(psec) 

14.3 

12.4 

14.1 

16.8 

13.7 

12.4 

Std  Dev 
(psec) 

15.1 

17.6 

16.3 

16.4 

15.4 

9.2 

This  testing  indicates  that  the  HAGR  units  can  provide  carrier  phase  observations  consistent  with 
a  time  transfer  performance  of  16  ps  1- sigma,  post- calibration.  The  testing  performed  using  the 
HAGR  highlighted  the  benefit  of  a  highly  stable  front-end  and  also  identified  key  requirements 
for  the  LO  generation  which  are  being  designed  into  our  core  product.  Testing  on  these  units  is 
continuing  to  show  their  phase  stability  from  turn- on  to  turn- on  and  also  repeating  these  tests 
over  temperature.  Testing  is  also  planned  using  a  dual- frequency  (L1/L2)  P(Y)  code  version  of 
the  HAGR. 


Time 


Time 


Figure  19  Unitl-Unit2  Time  Offset  Figure  20  Unitl-Unit2  Time  Offset 

(cycles)  Time  Set  1  (cycles)  Data  Set  4 


CONCLUSION 

In  conclusion,  the  test  data  taken  to  date  on  two  single- element  HAGR  units  indicates  that  the 
HAGR  is  capable  of  making  high  accuracy  pseudo- range  and  carrier  phase  observations  suitable 
for  precision  positioning  and  timing  applications.  Test  data  showed  that  the  pseudo- range 
accuracy  could  be  maintained  at  the  sub -meter  level  even  before  carrier  smoothing.  This 
significantly  improves  the  differential  accuracy  from  the  HAGR  and  also  speeds  the  time 
required  to  perform  carrier  cycle  ambiguity  resolution.  The  kinematic  test  data  showed  that  the 
relative  position  errors  were  within  4  mm  (horizontal)  and  6  mm  (vertical)  which  is  consistent 
with  a  carrier  phase  error  of  around  3mm  (1-sigma).  This  improved  accuracy  demonstrates  the 
advantage  of  the  HAGR  antenna  array  in  reducing  multipath  errors. 


The  carrier  phase  time  transfer  testing  performed  showed  that  the  carrier  phase  random  errors  (1- 
Hz)  were  maintained  at  around  16  picoseconds  (1-sigma).  When  these  errors  are  smoothed 
against  a  precision  clock,  the  time  transfer  error  could  be  expected  to  approach  the  tolerance  of 
the  HAGR  phase  calibration,  which  to  is  around  +/-  3  ps  in  these  tests. 

Further  improvements  in  pseudo- range  and  carrier  phase  multipath  reduction  can  be  achieved  by 
optimizing  the  HAGR  weights  to  adapt  to  detected  multipath  signals.  Since  the  HAGR  digital 


beam  steering  is  performed  under  software  control,  this  capability  could  be  added  to  the  current 
system  design.  Research  is  currently  being  conducted  on  the  performance  benefits  that  can  be 
achieved  through  an  adaptive  multipath  spatial  processing  algorithm. 

Based  on  these  results,  and  previous  testing  of  the  HAGR[4],  this  GPS  receiver  has  the  following 
advantages  for  precision  positioning  and  timing  applications. 

•  Highly  stable,  phase- coherent  front-end,  phase- locked  to  an  external  10  MHz  oscillator 

•  Increased  C/NO  to  the  satellite  observations  using  beam- steering 

•  High  accuracy  pseudo- range  and  carrier- phase  observations  for  rapid  carrier- cycle  ambiguity 

•  Multipath  minimization  on  both  pseudo-range  and  carrier-phase  from  the  digital  beam- 
steering 

•  L1/L2  P(Y)  code  HAGR  in  development 
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